Multi-gap superconductivity and Shubnikov-de Haas oscillations in single crystals of 

the layered boride OSB2 
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Single crystals of superconducting OsB2 [Tc = 2.10(5) K] have been grown using a Cu-B eutectic 
flux. We confirm that OsB2 crystallizes in the reported orthorhombic structure (space group Pmmn) 
at room temperature. Both the normal and superconducting state properties of the crystals are 
studied using various techniques. Heat capacity versus temperature C(T) measurements yield the 
normal state electronic specific heat coefficient 7 — 1.95(1) mj/mol and the Debye temperature 
©D = 539(2) K. The measured frequencies of Shubnikov-de Haas oscillations are in good agreement 
with those predicted by band structure calculations. Magnetic susceptibility x{T,H), electrical 
resistivity p(T) and C{T,H) measurements {H is the magnetic field) demonstrate that OSB2 is a 
bulk low-K [k(Tc) ~ 2(1)] Type-Il superconductor that is intermediate between the clean and dirty 
limits [C(r = 0)/£ = 0.97)] with a smaU upper critical magnetic field Hc2{T = 0) = 186(4) Oe. The 
penetration depth is A(r = 0) = 0.300 fim. An anomalous (not single-gap BCS) T dependence of A 
was fitted by a two-gap model with Ai(T = 0)/fcBrc = 1.9 and A2(T = 0)/kBTc = 1.25, respectively. 
The discontinuity in the heat capacity at Tc, AC/^Tc = 1.32, is smaller than the weak-coupling BCS 
value of 1.43, consistent with the two-gap nature of the superconductivity in OSB2. An anomalous 
increase in AC at Tc of unknown origin is found in finite H; e.g., AC/"/Tc ~ 2.5 for // ~ 25 Oe. 

PACS numbers: 74.10.-l-v, 74.25.Ha, 74.25.Bt, 74.70.Ad 



I. INTRODUCTION 



Although multigap superconductivity was first ad- 
dressed theoretically by Suhl et al. in 1959^ and the 
first experimental observation of the possible existence 
of two distinct superconducting gaps was made in 1980 
using tunneling measurements on Nb-doped SrTiOa,^ the 
subject of multigap superconductivity has only recently 
gained impetus after it was established that several un- 
usual superconducting properties of MgB2 could be ex- 
plained within a two-gap superconductivity scenarioj^ 
There are now several other candidates for multi-gap su- 
perconductivity like NbSe2 (Refjl), i^NisBaC {R = Lu, 
Y) (Ref.i), LuaFegSis (Refs.0,[l, and Sr2Ru04.^ 

In multigap superconductors distinct superconduct- 
ing gaps exist on different disconnected parts (sheets) 
of the Fermi surface (FS) although the interband pair- 
ing leads to a single critical temperature Tcjii^^ Most 
superconductors show multi-band conduction, but due 
to interband pairing the gap has the same magnitude 
on all bands. When interband pairing is weak then the 
gaps on different sheets of the FS can have significantly 
different magnitudes. This can lead to anomalous be- 
havior in the temperature-dependent heat capacity, up- 
per critical magnetic field iJc2, and penetration depth A 
measurement s j^i^ i ^ ^ d ^ 

The compound OsBa has a layered crystal structure 
qualitatively similar to that of MgBa, except that the 
B layers are corrugated in OsBa instead of flat as in 
MgBai^ The crystal strucure of OsBa is shown in Fig. [T] 
Figure [Ija) shows the crystal structure of OsBa viewed 
at a slight angle from the h axis. Figure [Hb) shows 
the structure projected on the ah plane. Along the 



c axis the boron layers lie between two planar transi- 
tion metal layers which are offset along the afe-plane. 
We have recently reportedi^ several anomalous behaviors 
for polycrystalline samples of the layered superconductor 
OsBa which has a superconducting transition tempera- 
ture Tc = 2.1 K.^^ These unusual behaviors include a 
reduced specific heat discontinuity at Tc in some samples 
and a magnetic field penetration depth versus temper- 
ature T dependence that was consistent with two-gap 
superconductivity. We also observed a positive curva- 
ture in the T dependence of the upper critical magnetic 
field Hc2- To gain further insights into these interesting 
behaviors, measurements on single crystals are needed. 

Herein we report the growth of OsBa single crys- 
tals, and structure, isothermal magnetization, dynamic 
and static magnetic susceptibility, specific heat, elec- 
trical resistivity, magnetic field penetration depth, and 
Shubnikov-de Haas (SdH) oscillation measurements on 
the crystals to characterize their superconducting and 
normal state properties. Following a description of the 
experimental details in Sec. [Hi the experimental results 
are given in Sec. IHIl A summary of the results and our 
conclusions are given in Sec. IIVI including a list in Ta- 
ble IIVI summarizing the parameters characterizing the 
physical properties that we obtained. 



II. EXPERIMENTAL DETAILS 

Single crystals of OsBa were grown with a high tem- 
perature solution growth method using Cu-B as the flux. 
First, a polycrystalline sample of OsBa was prepared 
by arc-melting Os powder (99.95%, Alfa Aesar) and B 
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FIG. 1: (Color online) The crystal structure of OsB2 viewed 
at a slight angle from the b axis. The Os atoms are shown 
as large (red) spheres while the B atoms are shown as the 
small (blue) spheres. A single unit cell (shown) contains two 
formula units, (b) Projection of the OsB2 structure onto the 
ab plane. 



chunks (99.5%, Alfa Aesar) taken in stoichiometric ratio. 
A Cu-B binary alloy was then prepared at the eutec- 
tic composition Cuo.87Bo.13 by arc- melting. For crystal 
growth the arc-melted OSB2 sample (~ 0.5 g) was placed 
in a 2 mL AI2O3 crucible. About 5 g Cu-B flux was 
placed on top of the OSB2 ingot. The crucible with a 
lid was placed in a vertical tube furnace which was then 
evacuated and purged with high purity Ar gas repeat- 
edly (« 10 times) after which the growth was started in 
a flow (« 60 cc/min) of Ar. The furnace was heated to 
800 °C in 30 min, then heated to 1450 °C in 6 hrs and 
held at this temperature for 6 hrs. The furnace was then 
cooled to 1020 °C at a rate of 2 °C/hr and then rapidly 
cooled to room temperature. Well-formed crystals with 
flat facets were obtained after the Cu-B flux had been dis- 
solved in dilute nitric (HNO3) acid. A scanning electron 
micrograph of a typical crystal is shown in Fig. [H^a). 

Some single crystals were crushed for powder X-ray 
diffraction (XRD) measurements. The XRD patterns 
were obtained using a Rigaku Geigerflex diffractometer 
with Cu Ka radiation, in the 29 range from 10 to 90° 
with a 0.02° step size. Intensity data were accumulated 
for 5 s per step. 

For single crystal structure determination, a well- 
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FIG. 2: (Color online) (a) Scanning electron microscope im- 
age of a typical OSB2 crystal, (b) Rietweld refinement of the 
powder X-ray diffraction data of crushed OSB2 crystals. The 
open black circles represent the observed X-ray pattern, the 
solid red line represents the fitted pattern, the dotted blue 
line represents the difference between the observed and cal- 
culated intensities and the vertical black bars represent the 
peak positions. 



shaped crystal (0.27 x 0.18 x 0.15 mm"^) was selected. 
The data collection for the crystal was performed us- 
ing a Bruker Apex II instrument with Cu Ka radiation 
at r = 100 K and was solved with latest version of the 
Apex software package which is reliable for a combination 
of numerical and multi-scan absorption correction. The 
initial cell constants were obtained from three series of uj 
scans at different starting angles. Each series consisted 
of 30 frames collected at intervals of 0.3° in a 10° range 
about uj with the exposure time of 5 s per frame. The 
obtained reflections were successfully indexed by an auto- 
mated indexing routine built in the Apex program. The 
final cell constants were calculated from a set of strong re- 
flections from the actual data collection. The data were 
collected using the full sphere routine by collecting 20 
sets of frames with 1 degree scans in w with an expo- 
sure time of 5 s per frame. This data set was corrected 
for Lorentz and polarization effects. The absorption cor- 
rection was a combination of a numerical one based on 
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TABLE I: Crystal data and structure refinement of OsB2. 
Here Rl = Ell^ob^l - li^caidl/El^'obsl and wR2 = (E[ 
w([J^obsP - |FealcP)']/E[«'(l^^obs^'l)'/^ where Fob, is the 
observed structure factor and -Fcaic is the calculated structure 
factor. 



Temperature 


100(2) K 


Crystal system, space group 


Orthorhombic, Pmmn 


Unit cell parameters 


a = 4.6729(3) A, 




b = 2.8702(2) A, 




c = 4.0792(3) A 


Unit cell volume 


54.711(7) A^ 


Z (formula units per unit cell) 


2 


Molar volume 


Id. 474(2) cm /mol 


Density (calculated) 


12.858 Mg/m^ 


Absorption coefficient 


212.327 mm-i 


F(OOO) 


172 


Data / restraints / parameters 


65 / 12 / 12 


Goodness-of-fit on 


1.015 


Final R indices [I > 2a{I)] 


Rl = 0.0294 




wR2 = 0.0736 


Extinction coefficient 


0.016(3) 



a face indexing and an additional correction based on 
fitting a function to the empirical transmission surface 
as sampled by multiple equivalent measurementsi^ using 
the Apex softwareJi 

The temperature dependences of the dc magnetic sus- 
ceptibility and isothermal magnetization were measured 
using a commercial Superconducting Quantum Interfer- 
ence Device (SQUID) magnetometer (MPMS5, Quantum 
Design). The resistivity and heat capacity were measured 
using a commercial Physical Property Measurement Sys- 
tem (PPMS, Quantum Design). The resistivity was mea- 
sured using a four-probe technique with a current of 5 mA 
along the b axis. The dynamic susceptibility was mea- 
sured between 0.5 K and 2.6 K using a 10 MHz tunnel- 
diode driven oscillator (TDO) circuit with a volume sus- 
ceptibility sensitivity Ax ~ 10"^!^^ The details of the 
measurement and the extraction of magnetic susceptibil- 
ity and penetration depth from TDO measurements have 
been described in our previous workJ^ 



III. RESULTS 

A. Crystal Structure of OsBa 

The powder XRD pattern for crushed single crystals 
of OsB2 is shown in Fig. [2{b). All the lines in the X- 
ray pattern could be indexed to the known orthorhombic 
Pmmn (No. 59) structure and a Rietveld refinemen1>i^ 
of the X-ray pattern, shown in Fig. [Hb), gave the lat- 
tice parameters a = 4.6855(6) A, & = 2.8730(3) A and 
c = 4.0778(4) A. These values are in very good agreement 



TABLE II: Atomic coordinates and anisotropic displacement 
parameters for OsB2 in space group Pmmn (second setting). 
The Os atoms occupy Wyckoff 2a (i,i,2) positions and 
the B atoms occupy 4/ (a;, i,2:) positions. Ui\, U22 and 
Us3 are the anisotropic thermal parameters in units of A^ 
defined within the thermal parameter of the intensity as 

^-2tv^ (h'^ a^Uii+k^b'^U22+l^':^Uaa) 



atom 


X 


y 


z 


Uii 


U22 


U33 


Os 


1/4 


1/4 


0.1527(2) 


0.003(1) 


0.004(1) 


0.003(1) 


B 


0.049(6) 


1/4 


0.359(4) 


0.004(8) 


0.001(8) 


0.006(8) 
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FIG. 3: (Color online) Electrical resistivity p for a single crys- 
tal of OsB2 versus temperature T with current 7 = 5 mA 
along the b axis and with applied magnetic field H — Oe. 
The inset shows the p{T) data between T — 1.7 and 2.6 K 
measured in various applied magnetic fields H as indicated. 



with our previously reported values [a — 4.6851(6) A, 
& = 2.8734(4) A, and c = 4.0771(5) A] for a polycrys- 
talline sample.^'' 

Single crystal XRD data were obtained at T = 100 K. 
The systematic extinctions of peaks in the XRD data 
were consistent with the space group Pmmn^ in agree- 
ment with earlier reports from single crystal and pow- 
der XRD measurements on OsB2.^° The positions of the 
atoms were found by direct methods and were refined 
in full-matrix anisotropic approximation. Some parame- 
ters obtained from the single crystal structure refinement 
are given in Table |T] and the final atomic positions and 
anisotropic thermal parameters are given in Table |lTl 
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B. Electrical Resistivity 

The electrical resistivity p versus T of a single crystal 
of OsB2 from 1.75 K to 300 K measured in zero applied 
magnetic field H and with a current 1 — 5 mA applied 
along the b axis, is shown in Fig. [3l The p{T) shows 
metallic behavior with an approximately linear decrease 
in resistivity on cooling from room temperature to 50 K. 
This behavior is similar to that observed earlier for a 
polycrystalline sample.^* At low temperatures p becomes 
only weakly temperature dependent and reaches a resid- 
ual resistivity po = 1.55 /if2 cm just above 2.2 K as seen 
in the inset of Fig. [S] The large residual resistivity ratio 
RRR = p{300 K)/po = 22 indicates a well crystallized 
sample. 

The inset of Fig. [3] shows the low T data measured in 
various H. The p in H = Q drops abruptly below 2.20 K 
and reaches zero by 2.14 K, as highlighted in the inset of 
Fig- El This superconducting transition was observed ear- 
lier by us for a polycrystalline sample;^ consistent with 
the original report in 1975 of superconductivity in OsB2 
by Vandenberg et al.^^ As expected, the superconduct- 
ing transition shifts to lower T with increasing H. These 
data were used to determine the upper critical magnetic 
field Hc2{T) which will be discussed later. In particular, 
for each applied field H , this H is taken to be Hc2 for 
the temperature at which the resistance drops to zero. 



C. Isothermal Magnetization and Magnetic 
Susceptibility 

1. Normal State 

The isothermal magnetization M versus H applied 
along the b axis is shown in Fig. IHa) at various T. At 
high T > 100 K, Af is diamagnetic and proportional to 
H with a slope that is almost constant between T = 100 
and T — 300 K. At lower T, M initially increases with 
H towards a positive value before showing saturation at 
a field of about 500 Oe. For higher H , M turns over 
and becomes diamagnetic. For > 1 T, M is linear in 
H with a nonzero y-intercept. A similar M{H) behavior 
is also observed in measurements with H applied along 
the a and c axes and is consistent with the presence of 
a small amount of paramagnetic and/or ferromagnetic 
impurities in the sample. Contributions from paramag- 
netic impurities are also observed at low temperatures 
in our normal state magnetic susceptibility x measure- 
ments via a Curie- Weiss-like upturn in Fig.[5]below. The 
M{H) data indicate that the impurity contribution sat- 
urates at high H. We extract the saturation magneti- 
zation A'Is{T) by fitting the M{H,T) by the expression 
M{H,T) = Ms{T) + XintH, where Xint is the intrinsic 
susceptibility. The Ms{T) data so obtained are shown in 
Fig.Hb). 

The normal state x = M/H versus T for OsB2, mea- 
sured between 1.8 K and 300 K with H — 3 T applied 
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FIG. 4: (Color online) (a) Isothermal magnetization M versus 
magnetic field H at various temperatures T with H applied 
along the b axis, (b) Saturation magnetization Ms versus T 
obtained by fitting the M{H) data above = 2 T at various 
T by the expression M{H,T) = MsiT) + x{T)H. 



along the a, 6, and c axes, is shown in Fig. [SJ The pow- 
der average susceptibility and the x(T) for a poly- 
crystalline sample from Ref. 14 are also shown in Fig. [5] 
The x(T) for single crystalline OsB2 is weakly temper- 
ature dependent between 50 K and 300 K. The upturn 
at low temperatures is most likely due to the presence of 
small amounts of paramagnetic impurities as mentioned 
above. The solid curve through the powder average x{T) 
data is a fit by the expression x(T) = Xint + j&g- The 
fit gave the values Xint — — 4.56(3) X 10-5 cmVmol, 
C = 1.43(1) X 10-4 cm^ K/mol, and 6 = -12.1(5) K. 
This value of C is equivalent to about 0.04 mol% of spin- 
1/2 impurities with a 17- factor g = 2. The large negative 
value of 9 is probably due, at least in part, to satura- 
tion of the paramagnetic impurities by the relatively high 
3 T field. The powder average value at 300 K is x(300 
K) = —4.50 X 10"^ cm'^/mol. A diamagnetic suscepti- 
bility for a transition metal compound is rare, but not 
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FIG. 5: (Color online) Magnetic susceptibility x versus tem- 
perature T for a single crystal of OsB2 with a magnetic field 
// = 3 T applied along the a, b, and c axes. The powder 
average x{T) the x{T) for ^ polycrystalline sample (from 
Ref. UM are also shown. The solid curve through the powder 
average x{T) data is a fit by the expression x{T) ~ Xint + j^- 
The intrinsic susceptibility obtained by correcting the pow- 
der averaged x{T) for the ferromagnetic impurity contribu- 
tion Ms{T), Xint(T) = x{T) - Ms{T)/H is shown as open 
diamond symbols. 
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FIG. 6: (Color online) Temperature T dependence of the 
zero-field-cooled (ZFC) and field-cooled (FC) dimensionless 
volume susceptibility Xv in terms of the superconducting vol- 
ume fraction (47rXv) of a single crystal of OsB2 with a mag- 
netic field H = 5 Oe applied along the (a) a axis, (b) b axis, 
and (c) c axis, (d) The T dependence of the ZFC super- 
conducting volume fraction 4nxv of a single crystal of OsB2 
measured in various magnetic fields H applied along the b 
axis. The construction used to determine Tc{H) is illustrated 
by the red straight line for H = 5 Oe. At low T, the 47rXv 
values are all more negative than —1 due to demagnetization 
effects. 



unprecedented 

Fig. [S] also shows the intrinsic susceptibility XintiT) ob- 
tained by correcting the powder averaged susceptibility 
for the presence of ferromagnetic and/or paramagnetic 
impurities as discussed above. Here we have assumed 
that Ms(T) is isotropic so that we can use the Ms{T) 
data obtained from M{H,T) data for the b axis to cor- 
rect the powder averaged susceptibility. 

The x along different crystallographic directions is 
anisotropic with the value (averaged between T = 100 K 
and 300 K) along the b axis (« -6.25 x lO'^ cm^/mol) 
being much smaller than the susceptibility along the a 
(« -3.75x10-5 cmVmol) or c (« -3.25x 10"^ cmVmol) 
axes which are quite similar. The similarity of the x val- 
ues along the a and c axes is surprising given the layered 
nature of the crystal structure which is built up of alter- 
nating Os and B layers in the a-b plane that are stacked 
along the c axis. However, theoretical Fermi surface 
(FS) calculations have shown that there are quasi-one- 
dimensional tubular structures running along the b axis, 
and the FSs along the a and c axes are quite similar i2i 

As described in Ref. [ij, one can estimate the para- 
magnetic Pauli spin susceptibility xp from the intrinsic 
susceptibility Xint according to 

Xint = Xorb +Xp, (1) 

where Xorb is the total orbital susceptibility, which in- 
cludes the diamagnetic core contribution, the paramag- 
netic Van Vleck contribution, and the Landau diamag- 
netic contribution from the conduction electrons. In 



Ref. [H, we estimated Xorb = -7.8 x 10-^ cm^/mol. 
However, the accuracy of this estimate for Xorb is un- 
known (see also below). Using this value of Xorb, our mea- 
sured Xint (300 K) = -4.50 X lO-'^ cmVmol and Eq. 
for our single crystalline OsB2 we obtain a powder aver- 
age Xp — 3.23 X 10^^ cm^/mol. 

From Xp one can estimate the density of states at 
the Fermi level N^ep) for both spin directions using the 
relatioEb^E 



XP = f^lNiep) = (^3.233 x 10"^ — j N{eF), (2) 

where fiB is the Bohr magneton and the equality on the 
far right-hand side is for iV(eF) in units of states/(eV 
f.u.) for both spin directions, where "f.u." means "for- 
mula unit." Taking the above average value of xp for 
OsB2, we get iV(eF) = 1.00(3) states/(eV f.u.) for both 
spin directions. This value is a factor of two larger than 
the value from our specific heat measurements below 
as well as from band structure calculations [N{eF) = 
0.55 states/(eV f.u.)],^'^ indicating that our estimate of 
the orbital susceptibility above is too negative. Using 
Eq. ^ and the band structure density of states value 
gives the revised estimate xp = 1-8 x 10"^ cm'^/mol. 
Then using the measured Xint 

(300 K) and Eq. O yields 
a revised powder averaged orbital susceptibility Xorh = 
-6.3 X 10-5 cm^/mol. 




H(Oe) 

FIG. 7: (Color online) Hysteresis loops at 1.7 K of the vol- 
ume magnetization My{H) normalized by l/47r, versus ap- 
plied magnetic field H applied along the (a) a axis, (b) b 
axis, and (c) c axis. The arrows next to the data indicate 
the direction of field ramping during the measurement, (d) 
Normalized magnetization M^[H) at various temperatures T 
versus H along the b axis. The construction used to deter- 
mine Hc2{T) is shown by the solid red line for T = 1.75 K. 
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FIG. 8: (Color online) Dynamic susceptibility x normalized 
by l/47r, versus temperature T at a frequency of 10 MHz 
with various applied magnetic fields H in units of Oe. The 
Atix data have been normalized to a minimum value of —1 at 
the lowest T. The construction used to determine Tc{H) is 
shown by the red line for H = 0. 



2. Superconducting State 

The temperature dependence of the anisotropic zero- 
field-cooled (ZFC) and field-cooled (FC) dimensionless 
dc volume magnetic susceptibility Xv of a single crys- 
tal of OsB2 measured from 1.7 to 2.5 K is plotted in 
Figs. [6lja), (b), and (c) in a field of 5 Oe parallal to the 
a, b and c axes, respectively, where Xv = M^/H and Mv 
is the volume magnetization. Complete diamagnetism 
in the absence of demagnetization effects corresponds to 
Xv = — l/47r, so the data have been normalized by l/iir. 
A sharp diamagnetic drop in the susceptibility along all 
three directions, below Tc = 2.05 K, signals the transi- 
tion into the superconducting state. The large Meissner 
fraction seen in the FC data for all three directions indi- 
cates weak magnetic flux pinning in the crystal. The data 
have not been corrected for the demagnetization factors 
Na {a = a, b, c) which give — 47rxv = j^j^ for the 
respective measured value. From the ZFC data at the 
lowest temperatures in Fig. [51 one obtains Na = 0.60, 
Nb = 0.60 and = 0.65, yielding Na + Nb + Ne = 1.85. 
This sum is greater than the value of unity expected for 
an ellipsoid of revolution. The reason for this discrepancy 
is not known. 

In Fig. ini^d) the temperature dependences of Xv mea- 
sured with various H applied along the b axis are shown. 
As expected the superconducting transition is suppressed 
to lower temperatures with increasing H. From these 
Xv{H) data the critical field Hc2{T) has been estimated 
using the construction in Fig. [6][d) , illustrated for H — 
5 Oe. The Hc2{T) has been determined by fitting a 
straight line to the data for a given field in the super- 
conducting state just below Tc and to the data in the 



normal state above Tc and taking the temperature at 
which these lines intersect as the Tc at that H. 

The hysteretic volume magnetization normalized 
by I/Att versus H loops measured at T = 1.7 K with H 
applied along the a, 6, and c axes are shown in Figs.[71[a), 
(b), and (c), respectively. There is a large reversible part 
in all the data recorded with increasing and decreas- 
ing H which again indicates very weak magnetic fiux pin- 
ning in the material. The M^{H) data recorded at var- 
ious fixed T with H applied along the b axis are shown 
in Fig.[7fd). Similar data (not shown) with H along the 
a and c axes were also recorded. The initial slope of 
the Mv{H) curves is larger than the value —1 expected 
for perfect diamagnetism, which indicates a nonzero de- 
magnetization factor, consistent with the 47rXv(T) data 
in Fig. ini From the My{H) curves in Fig. [TJd) we es- 
timated the critical field i?c2(T) from the construction 
illustrated in Fig. [ltd) for T = 1.75 K. 

The dynamic ac susceptibility x(^) measured between 
0.6 K and 2.5 K at a frequency of 10 MHz in various H is 
shown in Fig. [S) To determine 7Jc2 (T) from the data in 
Fig. Owe fitted a straight line to the data in the normal 
state and to the data just below Tc for a given applied 
magnetic field and took the value of the T at which these 
lines intersect as the Tc{H). This construction is shown 
in Fig. [5] for the data at H = 0. By inverting T^H) 
we obtain 7Jc2(T). The Hc2 has also been obtained in a 
similar way from the x(^) = M{T)/H SQUID magne- 
tometer data (not shown here) between 1.7 K and 2.4 K 
in various applied magnetic fields. 
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FIG. 9: Heat capacity C versus temperature T of a single 
crystal of OsB2 between 1.75 K and 48 K measured in zero 
magnetic field H. Inset (a) shows the C(r) data between 
1.75 K and 3 K measured with H — and H — 1 kOe applied 
along the c axis. Inset (b) shows the data for H = plot- 
ted as C{T)/T versus T between T = 1.75 and 2.5 K. The 
solid curve in inset (b) is a construction to estimate the heat 
capacity jump AC at Tc. 



The inset (b) shows the H = Oe heat capacity plot- 
ted as C{T)/T versus T between T = 1.75 and 2.5 K. 
The jump in the specific heat AC at the superconduct- 
ing transition Tc is usually normalized as AC/^Tc- From 
the construction shown as the solid curve through the 
data in Fig.lHlinset (b) we obtain AC/jTc = 1.32. This 
value is smaller than the weak-coupling BCS value 1.43/^ 
Considering the sharp anomaly observed at Tc it is un- 
likely that the smaller value of AC/jTc arises from a 
distribution of Tc's due to inhomogeneties in the sam- 
ple. We suggest that the small value of AC/"fTc arises 
from the multi-gap nature of the superconductivity as 
evidenced from our penetration depth measurements dis- 
cussed later. 

The electron-phonon coupling constant Acp can be esti- 
mated in the single-gap superconductivity approximation 
using McMillan's formula^ which relates the supercon- 
ducting transition temperature to Acp, the Debye tem- 
perature 8d, and the Coulumb repulsion constant ji* , 



9d 
1.45 



exp 



1.04(1 + Aep) 



Acp-/i*(l + 0.62Acp) 



(4) 



which can be inverted to give Acp in terms of Tc, Od and 
/i* as 



D. Heat Capacity 

Figure [9] shows the heat capacity C versus T data mea- 
sured on a single crystal of OsB2. The main panel shows 
the C(T) data measured in _ff = between T = 1.75 K 
and 48 K. The data below T = 10 K could be fitted by 
the expression C/T = 7 + /3T^ where the first term is 
the contribution from the conduction electrons and the 
second term is the contribution from the lattice. The 
fit (not shown) gave the values 7 — 1.95(1) mJ/mol 
and 13 = 0.0372(3) mJ/mol K^. From the value of /3 
one can estimate the Debye temperature 6*0 using the 



-vep 



1.04 + ^Mn(^) 
(1 - 0.62/1*) MtI^)- 1.04 



(5) 



expressioi 



.22 



From the value 0d — 539 K obtained above from heat 
capacity measurements, and using Tc = 2.1 K we get 
Acp = 0.41 and 0.50 for ji* = 0.10 and 0.15, respectively. 
These values of Aep are similar to those found for poly- 
crystalline samples and suggest that OSB2 is a moderate- 
coupling superconductor fi^ 

The density of states at the Fermi energy iV(eF) for 
both spin directions can be estimated from the values of 
7 and Acp using the relatiort^^ 



127r4i?n 



1/3 



(1.944 X 10*^) 



1/3 



/3J 



7 = 7o(l + Aep) 



(6) 



(3) where 



where R is the molar gas constant, n is the number of 
atoms per formula unit (n = 3 for OSB2), and the equal- 
ity on the far right-hand side is for 0d in K and /? in 
mJ/mol K^. We obtain Go = 539(2) K for OsBa. The 
values of 7, /?, and 0d obtained above are in very good 
agreement with the values we reported previously for an 
unannealed polycrystalline sampleJ^ 

The C(T) data below T = 3 K measured in zero and 
1 kOe applied field are shown in inset (a) of Fig. [5] A 
sharp step-like anomaly at T = 2.1 K is observed in the 
H ~ {) Oe data and confirms the bulk nature of the su- 
perconductivity in single crystal OSB2. The anomaly is 
suppressed to below T = 1.75 K in a field oi H ~ 1 kOe. 



70 = ^iV(eF) = 2.359 iV(eF), (7) 

A;b is Boltzmann's constant and the equality on the right- 
hand side of Eq. ([7]) is for 70 in mj/niol and iV(eF) 
in states/(eV f.u.) for both spin directions. Using the 
above 7 = 1.95 mJ/(mol K^), we find 7V(eF) = 0.59 
and 0.55 states/(eV f.u.) for the above Acp = 0.41 
and 0.50, respectively. These values are in excellent 
agreement with the value from band structure calcu- 
lations [iV(eF) = 0.55 states/(eV f.u.) for both spin 
directions] 1^ This agreement indicates that OSB2 is a 
weakly correlated electron system, consistent with the 
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-□- 15 Oe 
-A- 30 Oe 



6 7 

FIG. 10: (Color online) (a) Heat capacity C divided by tem- 
perature T versus T of a single crystal of OsB2 in vari- 
ous applied magnetic fields Ha.. A remanent magnetic field 
H m 25 Oe was present over the nominal value of the applied 
field Ha listed in the figure legend. The true H — Oe data 
from Fig. |9] inset (b) are also included for comparison, (b) 
C/T versus measured in various true magnetic fields H 
for a different OsB2 crystal. 



observed diamagnetic susceptibility in Fig. above. 
Henceforth we will take the bare density of states to be 

iV(eF) = 0.55 states/(eV f.u.) (8) 

for both spin directions, which corresponds to fi* — 0.15 
and Aop — 0.50. 

To obtain the critical magnetic field versus tempera- 
ture we have measured C(T) in various H. Figure [TUT a) 
shows the C{T) /T versus T data between T = 0.3 K and 
2.5 K, measured in various applied magnetic fields iJa- 
A remanent field of about 25 Oe was present in addition 
to the applied magnetic field Ha. in these measurements. 
For comparison the true H — Oe data from Fig. [S] in- 
set (b) are also shown. The superconducting transition 
seen as an abrupt jump in C{T)/T is suppressed to lower 
T with increasing as expected. However, the mag- 
nitude AC/jTc of the anomaly at Tc is initially larger 
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FIG. 11: (Color online) Upper critical magnetic field Hc2 ver- 
sus temperature T extracted from different types of measure- 
ments, as indicated. The straight line is a linear fit to the 
data near Tc- The dashed curve is a fit by the expression 
Hc2{T) = //c2(0)[l - i^)"]- The circle at T = K labeled 
"WHH" is the estimate of Hc2{T = 0) using the WHH for- 
mula in the clean limit (see text). 



than that observed in zero magnetic field. As shown by 
the arrows in Fig. rTOT a). AC/7TC increases from 1.32 
for _ff = Oe to 2.45 for _ffa — Oe (which is close to 
H — 25 Oe) suggesting a divergent nature of C at Tc in 
an applied magnetic field. 

The superconducting anomaly moves to lower T with 
increasing H and is not observed above Ha = 100 Oe 
{H sa 125 Oe). The step in aU data at about T = 1.15 K 
arises from a problem in the measurement and is not 
intrinsic to the sample. 

To further study the enhanced AC anomaly in low 
fields we measured C (T) for another single crystal in var- 
ious (true) H. The data are plotted as C{T)/T versus 
in Fig. [TUf b). We again observe that the anomaly at the 
superconducting transition becomes first order-like in a 
finite field showing that this feature is intrinsic to single 
crystalline OSB2. 

This AC{H) behavior is similar to that recently ob- 
served for r2Ga9 (r = Rh and Ir)^i2£ and for sin- 
gle crystals of ZrBi2,^'' where it was suggested that the 
Type-I superconductivity in these materials led to the su- 
perconducting transition in a finite magnetic field to be 
first order-like, resulting in a divergent AC{H) at Tc. A 
similar divergent AC{H) at Tc was observed 75 years ago 
for the Type-I superconductor thalliumi^ The C (T, H) 
behavior observed for single crystal OSB2 in Figs. [TUTa) 
and (b) is similar to that observed for the materials men- 
tioned above and might suggest that OSB2 is a Type-I su- 
perconductor. However, our estimates of the Ginzburg- 
Landau parameter k = 2(1) below indicate that OSB2 is 
a small-K Type-II superconductor. The unusual features 
in the C (T, H) for OSB2 are therefore not understood at 
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present but might be related to the multi-gap nature of 
the superconductivity. 



E. Upper Critical Magnetic Field Hc-iiT) 

The Hc2 (T) data obtained from all of the above mea- 
surements are plotted in Fig. 1111 In the temperature 
range of the SQUID magnetometer measurements (1.7 K 
to 2.4 K) all the data match well and the temperature 
dependence of Hc2 is linear (solid curve extrapolated to 
T = K in Fig. ini) with the slope ^ = -125 Oe/K. 
This linear slope can be used to get an estimate of the 
T = K upper critical field using the WHH formula for 
the clean limit i/c2(0) = -0.693 rc(^|Tj.— Using the 
above value of = -125 Oe/K and = 2.10 K we 

get ffc2(0) = 182 Oe. 

The Hc2{T) data at the lower temperatures T ^ 
1 K in Fig. [TT] show a deviation from linearity with 
a negative curvature. To obtain another estimate of 
Hc2{0), the Hc2{T) data in the whole T range were fit- 
ted by the empirical power law expression Hc2{T) = 



i?c2(0) 



1 



with -ffc2(0) and a as fitting param- 
eters and with fixed Tc — 2.15 K. The fit shown as 
the dashed curve in Fig. [TT] gave the values -ffc2(0) = 
188(2) Oe and a = 1.55(3). This estimate of i?c2(0) is 
close to the value of 182 Oe obtained above using the 
WHH formula. These two fits together yield our final 
value -ffc2(0) = 186(4) Oe. From Fig. [TT] it can also be 
seen that there is negligible anisotropy in the measured 
H^2{T) from T = 1.7 to 2.4 K. 

For a Type-II superconductor near Tc, the supercon- 
ducting coherence length ^ can be estimated from the 
measured Hc2 using the Ginzburg-Landau relation^ 



where (/>o = hc/2e = 2.068 x lO"'^ G cm^ is the flux 
quantum. We obtain an estimate of ^ using instead the 
zero-temperature value Hc2{T = 0) = 186(4) Oe arrived 
at above to obtain i{T = 0) = 0.133(2) jdia.. 



F. Superfluid Density 

The measured magnetic penetration depth A(T) in the 
superconducting state is related to the so-called London 
pentration depth Al (T) bj*^ 



A(T)«Al(T)Wi + 



(10) 



is the BCS coherence length, i is the quasiparticle mean 
free path and up is the Fermi velocity. Including the 
influence of £ gives the modified coherence length ^ as^ 
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1 



(12) 



The limit £^/£ is called the clean limit and the oppo- 
site limit the dirty limit. 

The superfluid density Ps{T) is related to A(T) by^ 



Ps(T) 



47re2A2(T)' 



(13) 



where m* is the effective mass of the individual quasipar- 
ticles, c is the speed of light in vacuum, e is the elemen- 
tary charge and ps is the density of quasiparticles that 
have condensed into the superconducting state, not the 
density of Cooper pairs which is a factor of two smaller. 
The normalized ratio of PsiT) to yOs(O) is simply 



(14) 



We now estimate whether OsB2 is in the clean or dirty 
limit or somewhere in between, by estimating the ratio 
£,{0)/£. The value of ^(0) was derived in the preceding 
section. We will estimate the mean-free-path £ using the 
measured resistivity at low temperatures and the A^(eF) 
in Eq. (|8|). First, the conductivity a is written as2^ 



ne^T 



(15) 



where n is the conduction carrier density and r is the 
mean-free scattering time of the current carriers. We 
then express t = vy/£, and from Eq. (|T5t we get 



h m*vpa 



(16) 



where in SI units the flrst term on the right is h/e^ — 
4108 n. Next we write both vp and n in terms of the 
known N{e-p) and then substitute these expressions into 
Eq. (HID. 

The (average) Fermi velocity has not been reported 
from band calculations. Therefore we calculate both 
vf and n from N(e-p) by assuming a three-dimensional 
single-band model with a spherical Fermi surface, 
yielding^^ 



where 



2m2(m*/me)' 



;V{eF), 



(17) 



Co 



hvF 

7rA(0) 



(11) 
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y 



2™ 



Mcf) 



(18) 
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FIG. 12: (Color online) The normalized superfluid density 

Ps(r)/ps(0) versus reduced temperature T/Tc- The solid 
curve through the data is a fit by the two-gap 7 model for 

superconductivity. The curves pi and p2 are the individual 
contributions from the two gaps. 
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FIG. 13: (Color online) Temperature T dependence of the 
two gaps Ai (top black solid curve) and A2 (bottom red solid 
curve). Also shown by the blue dashed curve is the BCS 
prediction for a single gap. 



where 2?(eF) is the density of states at the Fermi energy 
in units of states/(erg cm^) for both spin directions. Sub- 
stituting Eqs. (|T7)) and p8)) into ((T6)) . and using p = l/ci, 
gives 



3(m*/me)^ 



2m^ V ^ / 



(19) 



P 



where mo is the free-electron mass. The expression con- 
verting I'(eF) in units of states/(erg cm'^) to the conven- 
tional units of states/(eV f.u.) for both spin directons 
appropriate to the above definition of N{ep) is 



V{eF) = N{eF) 



1 1 / 1 eV \ TVa 

eV f.u. J V 1-6022 x IO-12 ergj Vm ' 

(20) 

where A'a is Avogadro's number and Vm is the molar 
volume. Substituting Eq. (|20l) into (|19p and putting in 
the values of the constants gives 



£ = 2.372 X 10 



_14 (m*/TOo)^Vf 



M 



VF = 2.622 X 10*^ 



NjeF) 
(m*/TOc)^VM 



(21) 



(22) 



where i is in cm, vf is in cm/s, N(eF) is in states/ (eV f.u.) 
for both spin directions, Vm is in cm'^/mol and p is in 
f2 cm. 

Inserting Vm = 16.47 cm'^/mol from TableU m* /m^ — 
1 (see SecHirHlbelow), N{eF) = 0.55 states/(eV f.u.) for 
both spin directions from our heat capacity data above. 



and p = 1.55 x 10"^ VL cm at 2.25 K from Fig. [3] into 
Eq. dHI) gives I = 0.137 pm at 2.25 K. Then using 
^(0) = 0.133 pm from above gives C(0)/£(0) = 0.97. 
Therefore OSB2 is in neither the clean limit nor the 
dirty limit, but in between. Irrespective of this difficulty, 
we will assume the clean limit in order to be able to 
carry out calculations for comparison with our measured 
penetration depth data. From Eq. ([^ we also obtain 
VF = 8.75 X 10^ cm/s. 

The ps(r)/ps(0) versus T/T^ calculated from the A(T) 
data using Eq. p4)) is shown in Fig. [121 The dashed curve 
is the T dependence of ps expected for a BCS single-gap 
s-wave superconductor. It is clear that our Ps{T) shows 
marked deviations from the single-gap BCS curve. This 
is consistent with our previous observations for polycrys- 
talline samples. The solid curve through the data is a 
fit by a two-gap 7 model^Si From the fit we obtained 
A(T 0) = 0.300 pm. 

The partial superfluid densities pi{T) and P2iT) from 
the fit to the two-gap model are shown as solid curves in 
Fig. [T^l The T dependences of the two gaps are shown 
in Fig. [13] plotted as normalized gaps A/fceTc versus the 
reduced temperature T/T^. For comparison, the T de- 
pendence of a single s-wave BCS gap is shown as the 
dashed curve. The T = value of the two gaps are 
Ai(0) = l.SSfceTc and A2(0) = 1.25fcBrc, respectively. 
The ratios of these two gaps to the single BCS gap value 
Abcs(O) = 1.76fcBTc are Ai(0)/Abcs(0) = 1.07 and 
A2(0)/Abcs(0) = 0.71. The values of these two gaps 
agree by construction with the theorem that in a two- 
gap superconductor, one of the gaps will always be larger 
than the BCS gap, whereas the second will always be 
smaller. 1° This constraint is a built-in result of the self- 
consistent solution to the two-gap 7 model. 
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G. Additional Superconducting Parameters 

The zero-temperature thermodynamic critical field 
Hc(0) of a superconductor is related to the zero- 
temperature superconducting gap A(0) in a single-gap 
BCS model by the expression'^'^ 

H.{0? P(eF)A(0)^ 

where, as above, 'D{€y) is the density of states at 
the Fermi energy for both spin directions in units of 
states/ (erg cm^). We use this expression as an ap- 
proximation to our two-gap model to obtain a value 
of ffc(O). Using the density of states value N{e-p) — 
0.55 states/eV f.u. for both spin directions from the 
above heat capacity measurements and Eq. ((20)) gives 
'D{e-p) — 1-26 X lO^'* states/(erg cm^). Using the larger 
gap found from fitting the penetration depth data and 
A(0)/fcBrc = 1-88 which gives A(0) = 5.45 x 10-^^ erg, 
Eq. ([23]) yields i7c(0) = 153 Oe. We can now de- 
rive the Ginzburg-Landau parameter n using the above 
ffc2(0) = 186 Oe via^o 



c2 



= 0.86. 



(24) 



This value is marginally on the Type-II side of the value 
K = 1/^/2 « 0.707 separating Type-I from Type-II super- 
conductivity, thus justifying the above notation of the 
measured critical field as being teh upper critical field 
Hc2 instead of the thermodynamic critical field He- 

Another estimate of k can be obtained using the 
relation^ 



n{T) 



[1 + {T/T,f 



2i/2 



27riJc(0)A(0)2 
0o[l + (T/Te)2 



(25) 



where k(0), i?c(0), and A(0) are the T — values of the 
Ginzburg-Landau parameter, the thermodynamic critical 
field, and penetration depth respectively. With the value 
Hc{0) = 153 Oe obtained above and the value A(0) = 
0.300 /im, we get k(0) = 6.0 and k{Tc) = 3.0. 

Two more estimates of k(Tc) can be made using the 
relations^S 
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FIG. 14: (Color online) The change A/ = f{H) - /(O) in the 
tunnel diode oscillator (TDO) frequency / versus magnetic 
field H measured at T = 180 mK with H applied along the 
c axis. The inset shows an image of the crystal used for the 
measurement. The c axis is out of the plane of the image. 
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FIG. 15: (Color online) The oscillating part of the TDO fre- 
quency df — A/— smooth background, versus the reciprocal 
of the magnetic field 1/H for H applied along the c axis. 



k(Tc) = 0.96-^^, (clean limit) (26) 



k(Tc) = 0.715-y^. (dirty limit) (27) 

Using A(0) = 0.300 /m, ^(0) = 0.133 ^m, £{T > = 
0.137 ^m andEq. (Uni), one obtains k(Tc) = 2.17 and 0.78 
from Eqs. (pS)) and (|77|) . respectively. 

The above four estimates of k(Tc) are all greater than 
I/V2 and therefore all indicate that single crystalline 
OSB2 is a small-K Type-II superconductor with k{Tc) = 
2(1). 



H. Shubnikov-de Haas (SdH) Oscillations 

Shubnikov-de Haas (SdH) oscillations in p{H) were ob- 
served as oscillations in the skin depth, which in turn 
were obtained from the oscillation frequency shift versus 
H oia. tunnel diode oscillator (TDO) in which the sample 
is placed inside the inductor of the LC circuit. Oscilla- 
tions were observed for T = 0.12-3.3 K in magnetic fields 
up to 7? = 14 T. Since ifc2(0) = 186 Oe from SeclHEl 
such fields quench the superconductivity and the mea- 
surements are therefore in the normal state. The inset of 
Fig. [13] shows an image of the crystal and the direction 
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TABLE III: Experimental (expt) and theoretical^i (thy) 
Shubnikov-de Haas frequencies F for single crystal OsB2 with 
the magnetic field applied along the c axis and in the ah plane, 
where n is the harmonic number. Also included are values for 
the effective mass m* jm^ and the electron-phonon coupling 
constant Aep of the electrons is specific orbits. For the data 
with H II ah plane, the field was oriented about 15° from the 
h axis, and hence the measured frequencies are not in general 
equal to the theoretical values for H \ h. 
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of the applied field H\\c axis where the c axis points out 
of the plane of the figure. 

The shift A/ = /(iJ)-/(0) in the TDO frequency ver- 
sus H applied along the c axis measured at T = 0.18 K is 
shown versus H in Fig. [TH Small oscillations can be seen 
riding on a smooth background. This _ff-dependent back- 
ground is due to the tunnel diode circuit that is partially 
exposed to the appHed field. The oscillations are more 
clearly visible when a smooth background is subtracted 
from Af{H) using a non-oscillating pieccwise cubic her- 
mite interpolating polynomial algorithm in Matlab. Fig- 
ure [12] shows the resulting oscillating part of the TDO 
frequency shift df versus the inverse magnetic field 1/H 
where df is the frequncy shift after the background sub- 
traction. 

To get the frequencies of the SdH oscillations at each 
T, a power spectrum was obtained by taking a Fourier 
transformation of the oscillation data such as in Fig. [151 
The resulting power spectra obtained for the measure- 
ments at T = 0.12-3.3 K are shown in Fig.[Tni The data 
reveal two clear fundamental frequencies Fi — 2767 T 
and F2 = 5905 T and possibly a third F3 = 660 T as 
marked in the plot in Fig. [161 although the intensity of 
the line at F3 is much weaker than the intensities of the 
prominent sharp lines at -Fi and i^2- We also observe 
the second harmonics for Fi and F2 but none for i^3, as 




Frequency (10^ Tesia) 



FIG. 16: (Color online) Power spectra of Shubnikov-de 
Haas oscillations obtained at the indicated temperatures from 
Fourier transformation of the quantum SdH oscillation data 
such as in Fig. [15] at the indicated temperatures. The spec- 
tra are shifted vertically by arbitrary amounts for clarity of 
presentation. 

shown in Fig. [16] and listed in Table ]lll] 

Recent first principles calculations of the Fermi surface 
(FS) showed three bands at the Fermi level, consisting of 
two nested deformed ellipsoidal surfaces (first and second 
bands) and a a corrugated tubular surface (third band) 
along the b axisj^ For a magnetic field applied along the 
c axis, the two closed electronic orbits which give rise 
to SdH oscillations are the cross-sectional areas of the 
two deformed ellipsoids normal to the applied field. The 
theoretically predicted frequencies of oscillations are 3023 
and 5902 Tj2i These values are in reasonable agreement 
with the experimentally observed frequencies Fi = 2767 
and F2 = 5905 T of the SdH osciUations for OsBz for 
measurements with H\\c axis. Since the frequencies of 
the SdH oscillations are inversely proportional to the area 
of the respective electronic orbits, we can assign the Fi 
oscillations as coming from the smaller inner ellipsoid 
while the F2 oscillations can be assigned to the outer 
ellipsoid. For the third tubular band, there are no closed 
orbits for H\ \c axis. The origin of the third frequency 
in Fig.[Tn]is therefore not understood at present. 

The T dependences of the amplitudes A{T) of oscil- 
lation in Fig. [in] can be used to estimate the respective 
effective masses for the bands responsible for the oscil- 
lations. The normalized A{T) is given by the Lifshitz- 
Kosevich formul a'^^'^^ 

and B H is the magnetic induction. The dimension- 
less variable X is proportional to the product of the ef- 



13 




15 



1.0 



I 0.8 



S 0.6 



■S 0.4 



I 0.2 
< 

0.0 



(b) 



m = 0.81 m. 



12 3 4 

T(K) 

FIG. 17: (Color online) (a) and (b) Fits to the amplitude ver- 
sus T data for the peaks Fi and F2 respectively, by Eq. (|28ll . 
The values of the respective effective masses m* for the two 
Fermi surface extrema are indicated in the figures. 



fective mass and the ratio of the thermal to magnetic 
energies. Therefore, by fitting the T-dependent ampli- 
tudes of the peaks in the power spectra in Fig. [16] by 
Eq. (pSj) one can obtain the m* values for the Fermi sur- 
face electrons responsible for the respective oscillations. 
FiguresfTWa) and (b) show A{T) for the Fi and F2 peaks, 
respectively. The fits by Eq. are also shown as solid 
curves through the data. We obtain rn{ — 1.03(4)TOe for 
Fi and rn^ — 0.81(3) mc for F2. The band masses pre- 
dicted by theory for these two orbits are mbi = 1.05mo 
and mb2 = O.SOtoc, respectivelyj^ The electron-phonon 
coupling constant Aop can be estimated by using the ex- 
pression m* = (1 -I- Aep)mb. Using the above values of 
TO* and TOb we obtain Aop ~ for Fi and Aop — 0.63(5) 
for F2. 

Figure [TH] shows the power spectra for quantum oscil- 
lations measurements performed with H applied perpen- 
dicular to the c axis. Due to the shape of our crystal we 
were only able to perform measurements with H tilted 
about 15 degrees away from the h axis as shown in the 
inset of Fig. 1181 The power spectra reveal nine different 
frequencies which are labeled in Fig. [TS]and their values 
are given in Table IIIII 

For H\\b axis, theory predicts at least six frequencies 
for quantum oscillations, two coming from the deformed 
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FIG. 18: (Color online) Power spectra of Shubnikov-de Haas 
oscillations obtained at temperatures T = 180 mK for H 
about 15° away from the b-cods in the a6-plane (see inset), 
obtained from Fourier transformation of quantum oscillation 
data (not shown) similar to those in Fig. 1151 The indices i for 
the quantum SdH oscillation frequencies Fi are as indicated. 
These frequencies are listed in Table IIIII 



ellipsoidal Fermi surfaces, and four from different closed 
orbits on the tubular Fermi surface as shown in Fig. 3 
in Ref. 21. For the two deformed ellipsoidal FS sheets, 
it was predicted that the frequencies of oscillation would 
vary linearly with angle as the field H is moved from H \ \ 
b axis toward 1 1 c axis<^ Thus, for the first deformed el- 
lipsoid the frequency of oscillation should change at a rate 
of about —25 T/degree and for the second deformed el- 
lipsoid it should change at a rate of about —3.8 T/degree 
as one moves away from || & axisi-^ The predicted val- 
ues of the frequencies of oscillation for H \\ b axis are 
5490 T and 6246 T, respectively, for the two deformed 
ellipsoids. If we use the above rates of change of the fre- 
quencies and the fact that we measured with H about 
15 degrees away from the b axis, then the expected fre- 
quencies are 5115 T and 6189 T, respectively, as listed in 
Table Hin These values are close to the observed frequen- 
cies Fq and Fg, respectively. Therefore, we can assign the 
frequencies Fq and Fg as coming from the two deformed 
ellipsoidal FS sheets. The assignment of the other fre- 
quencies in Table Hill is difficult since the angular depen- 
dence of frequencies arising from the tubular FS sheet is 
not known. We note that frequencies Fi and F2 in Fig.fTSl 
are much smaller and -Fg is almost a factor of 2 larger than 
any of the four frequencies predicted for the tubular FS 
sheet. We suggest that Fg is most likely a second har- 
monic of Fj. The remaining 4 frequncies F3 — 2983 T, 
Fi = 3812 T, F5 = 3957 T, and F7 = 5697 T can be com- 
pared to the theoretically predicted frequencies 2177 T, 
3265 T, 3888 T, and 5291 T^ The experimentally ob- 
served frequencies are similar to those predicted consid- 
ering the unknown angular dependence of the frequencies 
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arising from the tubular FS sheet. Thus we can tenta- 
tively assign these observed frequencies to the quasi-two- 
dimensional tubular FS sheet. 

The most intense frequencies in Fig. [THl F4, F5, Fq, 
and Fs , were used to estimate the effective masses by fit- 
ting the T dependences of these frequencies by Eq. (f28| . 
We obtain m| = 0.87mo, ml = 1.13mc, = 0.92mc, 
and mg = 0.96mo. The corresponding theoretically pre- 
dicted band masses are mb4 = 0.95mo, mb5 — 0.92mo, 
n^he = 0.62tocj and rribs — 0.45mo. Using the expres- 
sion m* = (1 + Aop)?Tib, we estimate electron-phonon 
interaction constants A4 « 0, A5 ~ 0.23, Ae « 0.48, 
and Ag ~ 1.13. It should be noted that the frequen- 
cies Fg and Fs arise from the deformed ellipsoidal FS 
sheets. Thus, we find that Acp is larger for the ellipsoidal 
FS sheets compared to the quasi-two-dimensional tubu- 
lar FS sheet. This suggests that the superconductivity in 
OSB2 is driven by the two deformed ellipsoidal FS sheets. 
The average value of Aop estimated above from McMil- 
lan's formula Eq. ([5]) was 0.4-0.5 which agrees with our 
inference that Acp is small on some FS sheets and is larger 
on others. 

The above experimental and theoretical SdH data are 
summarized in Table IIIIl 

IV. SUMMARY AND CONCLUSIONS 

We have grown high quality single crystals of OSB2 us- 
ing a novel Cu-B eutcctic flux. Measurements on these 
crystals confirm bulk superconductivity. The crystallo- 
graphic parameters of a single crystal are described above 
in Tables U and |TT] and Fermi surface properties in Ta- 
bic mil The various parameters describing other nor- 
mal and superconducting state properties are summa- 
rized here in Table HVl 

The heat capacity measurements show some unusual 
behaviors. The zero field anomaly at the superconduct- 
ing transition AC/jTc « 1.3 is smaller than the weak- 
coupling BCS value of 1.43. We suggest that this arises 
due to the two-gap nature of the superconductivity in 
OSB2. The occurrence of two superconducting gaps is 
supported by the anomalous temperature dependence of 
the penetration depth which could be fitted by the new 
7 model for multi-gap superconductors^ with the mag- 
nitudes of the two gaps being Ai(T = 0)/A:bT'c = 1.90 
and A2(r = 0)/A:bTc = 1.25, respectively. The zero- 
temperature upper critical field was determined to be 
Hc2{0) = 186 Oe. Four estimates of the Ginzburg- 
Landau parameter gave k{Tc) ~ 1-3 and thus indicate 
that OsB2 is a small-K Type-II superconductor. We ob- 
served an anomalous increase in the heat capacity jump 
at Tc measured in a finite magnetic field H. For example, 
at _ff « 25 Oe, AC/jTc « 2.5. This anomalous increase 
in AC/jTc was confirmed for two batches of crystals. 

The high quality of the crystals made it possible for 
us to study the anisotropy of the Fermi surface (FS) of 
OsB2 by measuring Shubnikov-de Haas quantum oscilla- 



TABLE IV: Parameters describing the physical properties of 
OsB2. The symbols are superconducting transition temper- 
ature Tc, electrical resistivity pb along the b axis, estimated 
mean-free path along the b axis just above Tc, £b{2.25 K), 
normal state magnetic susceptibility Xa {a — a, b, c) and 
powder average x measured in _ff = 3 T, linear specific heat 
coefficient 7, Debye temperature Qo, electron-phonon cou- 
pling constant Aep, bare band structure density of states de- 
rived from heat capacity measurements N{eF) for both spin 
directions, Fermi velocity wp, nearly isotropic upper critical 
magnetic field Hc2, superconducting coherence length ^, mag- 
netic penetration depth A, thermodynamic critical field He, 
and Ginzburg-Landau parameter k. 



Quantity 


value 




2.10(5) K 


Pb(2.25 K) 


1.55 nQ cm 


4(2.25 K) 


0.137 fim 


Xa(300 K) 


-3.9 X 10"^ cmVmol 


Xt{300 K) 


-6.3 X lO"'^ cmVmol 


Xc(300 K) 


-3.2 X lO"'^ cm^/mol 


X(300 K) 


-4.50 X 10"^ cmVmol 


7 


1.95(1) mJ/mol 


en 


539(2) K 




0.50 


iV(ep) 


0.55 states/(eV f.u.) 


v-p 


9.1 X 10'' cm/s 


i?c2(r = o) 


186(4) Oe 


ar = 0) 


0.133(2) nm 


A(T = 0) 


0.300 fim 


H4T = 0) 


153 Oe 




2(1) 



tions via contactless rf skin depth measurements. Some 
experimentally observed frequencies could be assigned to 
those predicted theoretically. The effective masses es- 
timated for the two deformed ellipsoidal FS sheets are 
larger than the predicted band masses and suggest a large 
electron-phonon coupling constant Acp ^ 0.5-1 for these 
FS sheets. A much smaller value of Acp was found for the 
third quasi- two-dimensional tubular FS sheet. These re- 
sults suggest that the superconductivity in OSB2 is driven 
by the two ellipsoidal FS sheets. This would also explain 
the negligible anisotropy in the measured upper critical 
fields between the three crystallographic directions. 
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